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ABSTRACT 

We analyse the stellar populations of 75 red-sequence dwarf galaxies in the Coma 
cluster, based on high signal-to-noise spectroscopy from the 6.5m MMT. The sample 
covers a luminosity range 3-4 magnitudes below M*, in the cluster core and in a field 
centred ldeg to the south-west. We find a strong dependence of the absorption line 
strengths with location in the cluster. Galaxies further from the cluster centre have 
stronger Balmer lines than inner-field galaxies of the same luminosity. The magnesium 
lines are weaker at large radius, while the iron lines are not correlated with radius. 
Converting the line strengths into estimates of stellar age, metallicity and abundance 
ratios, we find the gradients are driven by variations in age (> 6er significance) and 
. in the iron abundance Fe/H (~ 2.7<j significance). The light element (Mg, C, N, Ca) 

abundances are almost independent of radius. At radius of 0.4-1.3degree (~0.3-1.0x 
the virial radius), dwarf galaxies have ages ~3.8 Gyr on average, compared to ^6 Gyr 
near the cluster centre. The outer dwarfs are also ^50% more iron-enriched, at given 
luminosity. Our results confirm earlier indications that the ages of red-sequence galax- 
ies depend on location within clusters, and in Coma in particular. The exceptionally 
strong trends found here suggest that dwarf galaxies are especially susceptible to en- 
k> ; vironmental "quenching", and/or that the south-west part of Coma is a particularly 

. clear example of recent quenching in an infalling subcluster. 

c3 ' 
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1 INTRODUCTION 

Galaxies in rich clusters are subject to a range of environ- 
mental processes (stripping, suffocation, harrassment, tides) 
which are absent, or at least less efficient, in the field and in 
poor groups. The role of these effects in "quenching" star for- 
mation in cluster galaxies has been widely discussed based 
on instantaneous star-formation tracers such as Ha emis- 
sion. For example, the fraction of star-forming galaxies is su- 
pressed out to ~3 times the virial radius -R200 (e.g. Lewis et 

* Observations reported here were obtained at the MMT Ob- 
servatory, a joint facility of the University of Arizona and the 
Smithsonian Institution. 



al. 2002) . To trace the history of environment-driven quench- 
ing in today's galaxy clusters, high quality spectroscopy is 
needed to determine the final epoch of star-formation in 
now-passive galaxies. 

A number of works have claimed to observe significant 
variation in absorption line strengths of galaxies, as a func- 
tion of radius in clusters (e.g. Guzman et al. 1992; Carter 
et al. 2002; Smith et al. 2006). These trends have been in- 
terpreted as a radial gradient in average stellar age. Smith 
et al. note evidence that dwarfs may exhibit stronger radial 
dependences than the more massive galaxies, as expected if 
the efficiency of the quenching process depends on the depth 
of the potential well binding the gas reservoirs. 
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Figure 1. Sky distribution of the linestrength sample (black). 
The filled points are galaxies used in the fits of Section [3] Open 
symbols are either bright (M r < —19) comparison galaxies, have 
emission lines, or no measured age. Small points show other SDSS 
galaxies within the Hectospec survey region in red (confirmed 
members), yellow (confirmed background galaxies) or grey (no 
redshift information). The redshifts are from our Hectospec sur- 
vey and literature sources. Pale blue squares indicate the regions 
covered by HST imaging from the Coma ACS Treasury Survey 
(Carter et al. 2008). Crosshairs show the adopted cluster centre, 
midway between the two central cD galaxies. 



In this Letter, we provide further evidence for strong 
cluster-centric gradients in age for dwarf galaxies (median 
M r ps — 17.5) in Coma, using new spectroscopic data from 
the 6.5m MMT. A comprehensive analysis of the stellar pop- 
ulations will appear in a future paper (Smith et al., in prepa- 
ration). For luminosity and distance calculations we adopt 
a distance to Coma of 100 Mpc. 



2 OBSERVATIONS AND PARAMETER 
MEASUREMENTS 

High signal-to-noise spectra were obtained using the Hec- 
tospec fibre-fed spectrograph on the MMT, in February- 
April 2007. The instrument deploys 300 fibres over a 1 deg 
diameter field of view (corresponding to 1.75 Mpc at Coma); 
the fibre diameter is 1.5 arcsec (0.7 kpc). Two fields were ob- 
served, one centred on the cluster core, and an outer field to 
the south-west (Figure [TJ. The choice of the south-west field 
was motivated by available supporting data, and influenced 
by the work of Caldwell et al. (1993) who reported a high 
fraction of Balmer-strong galaxies in this region. 

Our observations were made in parallel with an exten- 
sive redshift survey of the cluster (Hornschemeier et al., 
in preparation). To study the stellar populations of dwarf 
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Figure 2. The g — r colour magnitude diagram for the galaxies 
in Figure [l] from SDSS photometry. For comparison, we note the 
position of the LF break, r*, and the magnitude range for "giant" 
early-type galaxies, having velocity dispersions cr>75kms — 1 . 
Symbols are as in Figure [T] 



galaxies, we observed 79 known cluster members with lumi- 
nosities 3 — 4 mag below M* , plus ten brighter galaxies for 
overlap with previous studies (e.g. Sanchez-Blazquez et al. 
2006). The target galaxies were selected to lie close to the 
red sequence of non-star- forming galaxies (Figure [2} . The 
270 line mm -1 grating was used, resulting in a wide wave- 
length coverage (3700-9000 A) at a spectral resolution of 
4.5 A, FWHM. The median total integration time for the 
faint galaxies was ~ 7 hours, yielding typical signal-to-noise 
ratio of ~40A~ 1 (at ~5000A). The brighter galaxies were 
observed for 0.7-2.0 hours. Relative flux calibration was im- 
posed using F stars with photometry from Sloan Digital Sky 
Survey (SDSS, Adelman- McCarthy et al. 2007), observed si- 
multaneously with the galaxies in each configuration. The 
data were reduced using hsred, an automated idl package 
based on the SDSS pipeline, provided by Richard Cool. 

Analysis of the one-dimensional spectra, including com- 
bination of multiple exposures and measurements of emis- 
sion and absorption lines, followed methods used by Smith, 
Lucey & Hudson (2007). Three galaxies with Ha in emis- 
sion (equivalent width > 20 A) were removed from the sam- 
ple. No other galaxies have emission above ~0.5A after re- 
moving the stellar continuum. Absorption-line indices were 
measured and corrected to the Lick resolution, but not trans- 
formed to the Lick "system" , since we will use models based 
on flux-calibrated spectral libraries. Comparisons with the 
index measurements of Sanchez-Blazquez et el. (2006) show 
excellent agreement for eight galaxies in common (e.g. HgF 
offset 0.1 ± 0.1 A with rms scatter 0.2 A). 

Line strength measurements can be transformed into 
estimates of stellar population age and element abundances, 
through comparison with population synthesis models. In 
this Letter, we analyse the absorption line measurements 
with reference to new stellar population models by Schiavon 
(2007). Advantages of this model set include its basis in a 
flux-calibrated stellar library, individually variable Mg, C, N 
and Ca abundances at fixed Fe/H, and a publicly-available 
code EZ-AGES for estimating parameters from measured in- 
dices (Graves & Schiavon 2008). 

The EZ-AGES code performs an iterative "sequential" 
grid inversion, finding abundance ratios which yield consis- 
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HdF +0.626+/- 0.143 (+4.4) 



CN2 -0.025 +/- 0.007 (-3.5) 



Ca4227 -0.136+/- 0.072 (-1.9) 
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Figure 3. Radial trends in the index residuals. The title bar for each index / gives the coefficient of projected radius R pro j in a bivariate 
fit of the form I = ao + a\M r + a2-R pro j. The significance of the radial trend is given in parentheses in units of the standard error. The 
red dotted lines show median residuals inside and outside 0.4 deg, while the dashed blue line shows an unweighted fit to the residuals. 



tent age and metallicity estimates across a range of index- 
index diagrams. A full description of the method is given by 
Graves & Schiavon (2008), while a summary can be found in 
Schiavon (2007). In our implementation, an initial estimate 
of age and Fe/H is made using Hbeta and the iron indices 
Fe5270 and Fe5335 assuming solar-scaled abundances. The 
Mg/Fe ratio is adjusted to yield the same metallicity and age 
from Hbeta and Mgb5177. Similarly, C/Fe is obtained using 
the Fe4668 index (which is heavily inflenced by the C abun- 
dance). With C/Fe in hand, the N/Fe abundance is adjusted 
for consistency with the measured CN2. Finally, Ca/Fe is 
obtained from Ca4227. The procedure is iterated, deriving 
a new age and metallicity estimate from the updated abun- 
dance pattern in each step. Once the final abundances are 
obtained, ages are estimated using the higher-order Balmer 
indices HdF and HgF, as well as using Hbeta. Errors are 
determined from Monte Carlo simulations. 



3 RADIAL TRENDS 

In this section, we examine the evidence for trends in our 
sample of Coma dwarfs, both at the level of observed indices 
and also using the derived ages and metallicities. Because 
the stellar populations of red galaxies follow strong scal- 
ing relations with the mass or luminosity (e.g. Smith et al. 



2007), we analyse the residuals from the index-luminosity 
(and age-luminosity relations, etc), rather than the index 
data themselves. Thus we are comparing the inner and outer 
galaxy population at fixed luminosity. (In fact, there is no 
correlation of luminosity with radius within our galaxy sam- 
ple, so similar results would be obtained without this con- 
trol.) We exclude the brighter comparison galaxies (with 
r < 16), and one object for which EZ-AGES failed to con- 
verge, leaving 75 galaxies in the fits. All galaxies are assigned 
equal weight in the fits, since the scatter is primarily intrin- 
sic and the measurement errors similar for all points. 

The residual index-radius relations are shown in Fig- 
ure [3] for the nine indices which enter into our analysis with 
EZ-AGES. Strong positive correlations are observed in all of 
the Balmer indices: Hbeta (5a), HgF (4a), and HdF (4a). 
Negative, but generally less significant, trends are recovered 
for the metal indices dominated by light elements: CN2 (3a), 
Mgb5177 (3a), Ca4227 (2a). Finally, the carbon-dominated 
index Fe4668 and iron-dominated Fe5270 and Fe5335 show 
no significant correlation. 

Figure [4] presents the equivalent relations for the de- 
rived age (estimated from each Balmer index, and the av- 
erage age), the metallicity Fe/H and the four light-element 
abundance ratios (Mg/Fe, C/Fe, N/Fe, Ca/Fe). We find a 
strong (5a) correlation towards younger ages at larger radii, 
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Figure 4. Radial trends in the stellar population parameters from EZ-AGES. Annotations as in Figure [3] 



whether measured with Hbeta, HgF or HdF. The size of the 
effect is consistent between age indicators. Using the aver- 
age of the three ages obtained for each galaxy, the age-radius 
correlation is significant at the 6.5a level, and its slope im- 
plies a factor of two change in age, per degree in radius. The 
median age for the galaxies beyond 0.4 deg from the cen- 
tre is 3.8 Gyr, corresponding to their star-formation being 
quenched at z<0.3 (Figured]). 

The iron abundance, Fe/H shows a significant (2.7a) 
correlation 0.20 ± 0.07 dex deg -1 , with higher metallicity at 
larger radius. (Note that the absence of a trend in Fe5270, 
for example, reflects the compensating effect of the age and 
Fe/H gradients.) The Mg/Fe abundance ratio is significantly 
smaller at large radius by 0.13 ± 0.04 dex deg -1 . Note that 
the magnesium abundance [Mg/H] = [Mg/Fe] + [Fe/H] is 
consistent with being independent of location in the clus- 
ter; thus the variation in [Mg/Fe] is driven by the greater 
iron abundance in the outer galaxies, not by lower magne- 
sium. The other abundance ratios are also smaller at larger 
radius, but with lower significance levels, ~2cr. The abun- 
dances N/H and Ca/H are, like Mg/H, consistent with no 
radial variation, while C/H may be marginally increasing 
with radius (2a). 

Magnesium is produced mainly by Type II supernovae, 
while iron is released primarily by Type la supernovae, 
which are delayed with respect to star formation. Therefore 
the smaller Mg/Fe (and larger Fe/H) at large radius sug- 



gest the outer galaxies had extended star-formation histories 
prior to their quenching. In this scenario, it is qualitatively 
expected that a trend to younger ages would be accompa- 
nied by a trend towards lower Mg/Fe. Quantitatively how- 
ever, galactic chemical evolution models are required to de- 
termine what star-formation and enrichment histories could 
reproduce the observed behaviour. 

Finally, we note there is a significant radial colour cor- 
relation for the sample, with a slope —0.05 ±0.01 mag deg -1 
in g— r (bluer colour at larger radius). A small shift, with the 
same sense, was previously noted by Terlevich, Caldwell & 
Bower (2001) and Carter et al. (2002). The galaxies beyond 
0.4 deg have higher recession velocities by 340±220kms -1 
than those in the core. Preliminary profile fits to SDSS im- 
ages do not reveal significant radial trends in morphology. 



4 DISCUSSION AND SUMMARY 

Several authors have claimed to observe significant variation 
in absorption line strengths as a function of radius in clus- 
ters, with much of this work focused on the Coma Cluster. 

Caldwell et al. (1993) were the first to note a large frac- 
tion of galaxies with strong Balmer absorption in the south- 
west part of Coma. Their sample probed a smaller radial ex- 
tent (out to ~50 arcmin) than our work, and was restricted 
to brighter galaxies (r<16). Their result has driven other 
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Figure 5. Quenching redshift, denned here as the redshift for 
which the look-back time is equal to the Balmer age, as a func- 
tion of projected distance from the cluster centre. Although the 
most recent burst will generally be superposed on an older popu- 
lation, the Balmer age is strongly weighted toward the youngest 
stars present. If the mass-fraction of the most recent burst is very 
small, the quenching redshifts plotted here are upper limits. The 
assumed cosmology has parameters (Qm ,&A, h) = (0.3,0.7,0.7). 

work in the south-west region, influencing the field place- 
ment of subsequent studies including the current one. 

Carter et al. (2002) observed Coma galaxies over a wide 
luminosity range, r ~ 13 — 19, in the cluster core and an 
outer field in the south-west. The fainter part of their sam- 
ple corresponds to typical luminosities in our MMT work, 
and the radial extent is comparable. They found a signifi- 
cant radial decrease in the Mg2 index at fixed luminosity, 
interpreting the result as a metallicity gradient, with lower 
metallicities in the cluster outskirts. However, their data for 
iron-dominated indices do not share the radial trend, while 
Hbeta is found to increase with radius. Taken together in 
fact, the results favour an interpretation in terms of age, 
rather than metallicity. Figure [6] demonstrates good agree- 
ment between Carter et al. (2002) and the present work, 
at the level of the index gradients. Our independent data, 
of much higher signal-to-noise, strongly confirm the earlier 
results for dwarfs in the south-west part of Coma. 

For giant galaxies, outside of Coma, the situation ap- 
pears quite different. Smith et al. (2006) analysed line- 
strengths for ~3000 galaxies in ~90 clusters of the NOAO 
Fundamental Plane Survey (NFPS). The NFPS is domi- 
nated by M* galaxies, i.e. with luminosities 15-40 times 
larger than our MMT sample. Smith et al. found stronger 
Balmer lines and weaker light-element features (Mgb5177, 
CN2, Fe4668) at larger distance from the cluster centres. 
The iron-dominated indices show no radial dependence. This 
pattern can be reproduced by variations in both the age 
and the a-element abundance ratio [a/Fe] (comparable to 
[Mg/Fe]), with no change in the overall metallicity [Z/H] 
(comparable to [Mg/H] here). 

The NFPS index gradients are much shallower than 
those found here or by Carter et al. (Figure [6} . The rel- 
ative pattern of gradients is also different, with stronger 
gradients in the light-element indices than in the Balmer 
lines. This is reflected in the gradients of derived param- 
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Figure 6. Comparison of our index gradient pattern (yellow) 
with previous work. For clarity, the results for each index have 
been normalised by its age response. Thus if age was the only 
variable changing with radius, all indices should give the same 
gradient. Yellow boxes show the la range from this work. For 
Carter et al. (2002, green) we compare their (Fe) with our Fe5335; 
their Mg2 gradient has been converted to an equivalent gradi- 
ent in Mgb5177. For Smith et al. (2006, red), we have converted 
their gradients in R200 to degrees assuming R200 ~ 80 arcmin for 
Coma. The bottom row shows the age gradients derived for the 
three studies, including a new estimate for Carter et al., based on 
their published index gradients. 



eters: the Hectospec age gradient (0.35 dexdeg -1 ) is seven 
times larger than than the NFPS trend (0.05 dexdeg -1 ), 
while the Mg/Fe trend (0.12 dexdeg -1 ) is only a factor of 
2.4 steeper than the NFPS trend in [a/Fe] (0.04dexdeg -1 ). 

In conclusion, we have confirmed earlier indications that 
the spectra of faint red-sequence galaxies in Coma are cor- 
related strongly with their distance from the cluster centre, 
at least toward the south-west. Our analysis, using the lat- 
est stellar population models, provides a clear interpretation 
of the observed trends as being primarily due to age. Since 
such strong radial trends are not seen in general for giant 
galaxies in clusters, a key question is whether the difference 
arises only from the different luminosity ranges, or whether 
the south-west region of Coma is unrepresentative. If the lat- 
ter, we may be observing the descendents of star-bursting 
galaxies in cluster-feeding filaments at z ~ 0.2 — 0.3 (e.g. 
Fadda et al. 2008). Future work must test whether dwarfs in 
other clusters, or in other parts of the Coma cluster exhibit 
similarly strong cluster-centric correlations. As one hint that 
the effect may not be confined to this region, Michielsen et 
al. (2008) have observed a sample of 18 dwarf ellipticals in 
Virgo, and find only young (~2 Gyr) objects beyond a ra- 
dius of ~1 Mpc, equivalent to ~0.6 deg at Coma. A definitive 
answer for Coma itself should result from an extension to 
our Hectospec programme, which will obtain data in more 
fields to provide a complete azimuthal coverage to a radius 
of 1.5 deg, or ~2.5Mpc. 
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